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1.0 INTRODUCTION 


Included in this report is a novel method of inverting Laplace transforms by 
using a new set of orthonormal functions. As a byproduct of the inversion it 
is seen how to approximate very complicated Laplace transforms by a transform 
with a series of simple poles along the left-half plane real axis. The inver- 
sion and approximation process is simple enough to be put on a programaolc hard 
calculator. 

2.0 INVERSION AND APPR OXIMATION 

Let f(a) be a Laplace transform and F(t) its exact inverse. N P(t) will be 
the approximate inverse, given by . 


*p(t) - A*] L ] (st ) + AgL-^C st) 4* •• • *• Aj^Lj^(st) 


0 ) 


where the Ly,(st) are the new orthonormal functions (described below and in the 
appendix). The A n values .are the Fourier coefficients and are given by 


" <s) ■ I r 


•’( t)L n (st )dt 


(2) 


3 is a free parameter chosen to produce the best approximation, as shown below. 
The integral square approximation error is given by 


E(s) 


^F(I) - F(t) y'dt = / F ( t ) *■ d t 


N 

V 

n=1 


A n (s) > 0. 


( 3 ) 


To minimise the integral square error, s is chosen such that 


N 

V 
*■■*»»*' 
nr 1 


A n (o) is maximum 


(M 


/ 
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The new orthonormal functions are shown below. 

t n = nai e- st + n a 2 e" 2st + n a 3 e ' 3st + ••• + n a n e " nst 
The values of n a^ are chosen such that 

CO 

f = 0 for n i m 

0 =1 for n - m 

The first 10 orthonormal functions are listed below. 

L i = /?se” i;t s > 0 

Lo * /M-re- st + 3e“ 2at ) 

L 3 = ,'T>s(3e“ st - 1.?e~ 2st + 10e-3st) 

L)i s ,'8s( _!i,»-st + TOe" 2 ^’ - 60e”' r ‘^ + 35e _4st ) 

L r - - /T0a(‘-o- 3t -60e“ 2st ♦ 210e"? st - 280e” ,|st + 126e”5st) 

l : /T2s(-be“ st + lObe”'-' 3 '' - 560«;“3st + I260e“' tat 
- + i|6,'_'e"' r! ) 


( r >) 


( 6 ) 


L 7 - /fITs(7c- st - I68e~'- St + 1260e~3 3t - il200e~ 4st + 6930e“5st 
- 55i»lie“ 6at + 17l6e- 7st ) 

Lg r /iFs(-8e“ st + 2‘32e“-' >3t - 25?0e-3»t + 1 1 ‘3<30e“ l4st - 27720e'5 st 
+ 360?6e“ 6st - 2'lo;'!ie-'i' st + 6»l35e~ 8st ) 

, i0 . /7Fi(<v-st . 3 ( ;,0o + »l620o”33t 27720e~ ,|at + 90(K!O C -'3st 

_ i68it>8e- fc:it + I02960e~ 8st + 2H310e“*; t ) 
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L 10 • v'25s(-10e-st + H95e-2at - 7920e*3st + .60060e-“»t _ ?52252e“5st 
+ 630630e- 6st - 960960e-? st + 875l60e- 8st ' - iJ37580e-9st 
* 92378c-10st) 

Figure 1 shows plots of the first four. L r> . 

The values of the n a< coefficients are given by 


| nai = <-1) ri+i /2 


sn- 


(n + i - 1 ) f 
i ! < 1 - 1 ) I (n - i)! 


(7) 


or 


n a j = (-1) n+i *'2sn 


i I ( i - D! j=1 


1 n , (n 2 - j 2 ) 


( 8 ) 


where n a-j = (-1) n+1 * / 2ns n 


(9) 


The recursion relationship for the n a j is given by 


n a l 


= (-1) ,1+ ^ •■7ns n 


n = 1, 2, , N 


( 10 ) 


n'H - 


n 2 - ( i - 1 ) 2 

i’( i l’)’’ 


i ~ 2 , 3 , • h 


(ID 


n = 2, 3, N 


The recursion relationship for the L n is given by 


2n - 1 
*n(n - 1) 


n = 2,3, •••, N 


( 12 ) 


n ( n - 1 ) 


cn - i 


n = 2 , 3 , ‘ • • , N 


(13) 
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L-j = /2s e" st 

= /s e“ s ^(6e“ at - ft) 

L n « U n ((e-«t - V n ♦ V n 

This is the relationship that should be used to compute the L n in a computer 
program. It is simple, fast, and accurate. 

i 

I - ■ 

Equation 2 'gave the Fourier coefficients in terms of F(t). In tt?rn:s of 
the Laplace transform, f(s), they are given by 


A n = N r 3if(is) j 

1 T 1 ‘ I 


Note that as n increases, so does the magnitude of the n a i» which has 
m oscillating sign. This can cause serious roundoff error problems in 
computing the A n . It is speculated that the maximum value of n = N be 
limited to approximately the number of significant decimal digits of accuracy 
used by a particular computer. One way to evaluate- this problem lor a particular 
e o mp u t e r 1 s> t q set' 1 


(1ft) 

05) 

_l)L n _i ~ L n _ 2 /U n _0 n = 3, *», N ( 16 ) 


f(s) 



Let s = 1 a ! \ d oo mpu te t h c A rJ . Tit eo re t toall y 

A i = 1 /V? ’ 

A., = 0 for n > 1 


‘" J Also see theorem 15 in the 



' «r ww r 




and 
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N 

^ r ^ 55 ^ ^ 

n= 1 


Due to roundoff error, the theoretical values will not be achieved for N 
large. 


Perhaps a better way of computing the A n (which may be slightly less affected 
by roundoff error*) is to use the algorithm shown below, which also computes C. 

C s 0 

DO c n s 1 , N 

Ap ~ f ( r. s ) 

IF (n.FQ, 1 ) GOTO b 


6 s 1 

DO a i r !, n - 1 


n n s 


* ( - n - 1> , • f 

> ■ n - i ; 11 


( n + 1 « i K n 


a 6 = ~6 


A n = r'?A )3 nA n 


C = C ♦ Ar 


c PRINT C 



^ .M W?? ?QW'**Ql ?K XfWpMft 


1 f #}**&* ?.'.• * ^7*t**- 
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Note that a should be ch os en su ch that C i 3 maxi m* j m * 

All the L n approach zero as t approaches infinity. Therefore, the 
approximations work well only when F(t) -► 0 as fc approaches infinity. This 
will be the case for stable system weighting functions - an important applica- 
tion. An example of what to do. when F(t) does not decay to zero is shown 
below. Let 


g(s) s 


1 - e- 2fi 




Apply the final value theorem. 


G<*} s Lin sg(s) s 2 
3*0 


3o instead of inverting g(s), invert 


f(s) = g(s) - ~ 


-i 


Now F ( t ) * 0 as t approaches infinity and G(t) - K(t) * 2. Thus 
z 2 + *\F ( t 1 

3« 0 EXAMP LES 

As the first example, let 


3 + " 

f(s) = r 

( S +• 1 ) t: <*■ IT'- 


( 18 ) 


The exact inverse is 

F(t) = e“^ cos (nt) 


v. 
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Figures 2 through 9 show the values of 


N 

c» y 

n= 1 


( 20 ) 


versus a for valuer, of U from l to Vi. The maximum value that C can 
obtain (neglecting roundoff errors) is 0.27300 since 


/ 


■(t) 2 dt = 0.27300/ 


(21) 


It is s 

ic<sr that each value of fi 

choice 

of 3 can greatly influent 

t 

Figure 

I 

10 shows plots of! F( t } , 

of 3 

was a - 2 . 2 . 1:; 1 'his case 


'A 5 s 0.33378 95910 


A p - 0.28719 57089 


A 3 - -0.21931 58987 


'.as its own optimum value of s, and the 
? the accuracy of the fit. 

(t)» p ' r 'd fjF(t). For N = J the optimum valu* 


3 F < t. ) s •3o^ c 3e* 2, ^ t + 1 i . 92 le~^ • /;t « 7.805e“6.6t 
liie approximate Laplace transform is thus seen to be 

3.3-5 • 11.Q21 


3 f(3) 


.V805 


s + 6.6 


For ft s 6 the optimum value of 3 ± 0,0 and' 

A * - .0.1 5 9 * 0 f 2 2 1 S A ^ v 0.3 6 f 0 7 f 1 7 f 


( 22 ) 


( 23 ) 
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A 3 S 0.22324 71254 A 4 * -0.14742 23756 

A 5 * -0.11047 63746 Aft * 0.11492 41046 

ftFU) 8 -1.9l6e“°*9t ♦ 43.627c"' 1 - 8t - 252. 178e*2.7t 

♦ 554.631c-3-6t . ‘il7.636e-4.5t * I74.488c~5.4t (;>n) 


From figure 8 it iu seen that N s 


fit. For ' 

:hta case 



Ai 

= 0. 1 ^9^0 

23073 

A 3 

= 0*31771 

101 

0 3 

As 

r -0. 16661 

71 

768 

A? 

= O.IOJCM 

315 

nS 

Ad 

~ -0,03666 

09 

018 


10 and s * 0.65 will give an excellent 

A"> s 0.31 134 50651 
Ah s 0.03684 84056 
Aft = -0.03590 26947 
Aft r -0.02596 24761 
\ 5 0 = 0.04125 66606 


and 


1 0 F ( t ) s -0.832e-°’6‘5t 59.853e-1.3t _ 1 195.825e~ 1 *95t 

♦ 10 136. 3?4e“^«6t _ u'4 250.C68e~3.25t * no 05o.528e~3.9t 

- 169 u40 . 63 ?.e~4 . 5ut + 1 42 . 52.6 • 1 34e“5 • *■ 6 

- 68 134.644e~5.85t . 13 742. 1?l e -6.5t (25) 


For the next exan.^le 


f(ij) s In 


{.:ti 


rhe exact inverse 1 a 



I 

! 


F(t) = -(e-t - e -2t) _ e -2t 
t 
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( 27 ) 


Figure 11 shows plots of C = 


t 


2 

A n versus s for N = 4, 8, and 12. 


■It is 


clear that for 
In this case 


s = 0.5, only four terms are needed to give an excellent fit. 

i 


4F(.t) = -0.0026le~°-5t 


Figure 11 shows plots of F(t) 
between F(t) and i,F(t). 


+ p- i7291e -t + 0.663l6e" 1 *5t _ 0.83353e~2t ( 28 ) 

I 

I 

I 

and nF(t). There is no visible difference 


Let 


g(s) = In 


(s + 2 
+ 1 


(29) 


Then from equations 26 and 28, g( s ) is approximated by 

. pfal . °- 00261 0.17292 0.66316 0.16647 

ll o ' 0 * - “ + — + , 

s 4- 0.5 S + 1 s 4 1,5 S 4- 2 (30) 

For s > 0, „g(s) is an excellent approximation of g(s), as seen below. 


9 



y >’<■>-< * •►.•jS'V; )~t vvj* -e. • 


M {!<;r>W!rr5^>'v' r? ; *iy -» *»r-T 


-.-• . TfJVV fV-. V* ft w* ^ r r F i ■;•■ « -f>' Vft - l/r**** y.^ 
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3 

g(s) 

2,8(3) 

0 

0.693D 

0.69304 

0.1 

.64663 

. 64660 

1 

. 1105*47 

. 11054? 

2 

. 28768 

. 28769 

5 

. 15415 

. 15415 

10 

o 

c— 

CO 

o 

. 0870 1 


Note 


G(t) = -0.00261 


o -0.5t + o.17292e -t + 0.663l6e _1 - 5t + 0. 166*176' ~ 2t (3D 


where 


G(t) = - (e“ l - e _>t ) 
t 


(32) 


Note G( 0 ) = 1 ami i|G(0) = 0.99994, 
For tlie final example 

f(s) = e" 

which has an exact inverse of 


(33) 


1 

F( t) = — - exp 

'VnO ^ 



(3'0 


Figure. 13 shews plots of 0 - \ A n versus s for values of N = 6, 10, and 1*1 . 


n= I 


Kor N •- 6 the optimum value of a ~ 0.8, and 


10 







ORIGINAL PAGE IS 
OF POOR QUALITY 





0 


] 






s 


Figure 4.- 


z a; 


f( 1 5 - 


e -t cos.(jrt). 








r " w i. / ^.w>,7r.r:-!^g w B ^y.wr;» 7 ^ T 'y^ ww^jy^HW^C My! 



















APPENDIX 


80rM.':0 


THE i, R FUNCTIONS AND PROPERTIES 

For brevity, the theorems) and lemmas presented here win be shown without 
proof. • 

Definition 1: 

The scalar product of f(») and g(t) will be defined by 


(f,g> 



f( t)g( t)dt 


Definition 
•'-fine L n Ut) by 


O) 



r l 


( 2 ) 


••’here 


r; M = (-1 v ' l V. 


Cn + i - n * 

I - ~ ~1 


\ :> ) 


A i t “mat o ly 




1 ) ! 


1=1 


• ‘O 


where 


„ v>: 


r: 1 - 


a 
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Lemma 1 : 
For n > 1 


" n a i (x - 1)(x - 2) • ■ • (x - (n - D) 

A iTTi ’ ' Sn (x + 1)(x + 2) • • • (x + (n - 1))(x + n) 
i=1 • 


( 6 ) 


Corollary A: 


/ 


L n (st)dt = (-1 ) n+1 j 


!• 2 


or 


/5ns 


(7) 


n a i 1 _ 1 

% — = (-i)n+VSsH - 
^ i n 

i=1 


(B) 


Corollary B: 

L n (0) = /2nJ 


(9) 


or 


n 



i= 1 


23 



Theorem 1 : 

The system of functions l, n (st) are orthononnal. 

(L n ,L m ) = 0 for n i m 
= 1 for n = m 

Definition ! t : 


The generating funetion g(z,t) is defined as 

— = g(1/z,t) 


g(z,t) = 1 


/ 


Tz 


\r 


*-st 


(1 - z) : 


Theorem 2 : 

Expansion of g(z,t.) into Maclaurin's series give 


P(a.t> = 


W 1 

\ st) 


< 1 


n=1 


/2ius 


BU.t) 




V 

nTl 


(- 1) n+1 

/2ns. 


1 

— L„(st) 


-n 


r p 


Theorem 3: 

T he d i f fe r e n c e eq u a t i on sa t .is f ' • d by L n ( s t ) 


That 


is 



■ IV *T|'- •* 






L n = <? 


2n - 1 ( 


/nCn^D (I 




3 t _ 


n ( n - 1 ) ( n - 1 ) ( n - 2 ) 


2n - 1 


2n 


' n " -'H , 

7T - j L „.1 


1 yTn - 1 )(n - 2 )^ ^ 

2 2n - 3 n "”‘ j 


Theorem ! l : 

The differential equation satisfied by L n is 

1 . 

|( e st - 1 ) L n + ae st L n .* s 2 n ? I, n = 0 
Also of interest is . 


L n * (-D 


n+ I 


d n-1 


(n - D ! d ( e “ st ) 


f e - nst (1 - e - st ) n_1 ] 

n-1 L 


iheorom 5 : 


e -nat . : „!( n - 1)! 


,v 


/T 


i=i 


( n + i ) J ( n - i ) 1 


■L] (at) 


Pof . init.ion b : 

i.Ot 


F ( t ) * ■ d t 


bo finite. 


Lot jy|F(t) bo an approximation of' F(t). 
do fined by 


The integral square error i 
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( 1 b) 


(16) 


( 17 ) 


(18) 
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E= J( N F(t) - F(t)) 2 dt 

0 


(19) 


Theorem 6: 

■ . - Rm in the integral square error sense 

The best approximation of F(t) u\ ™e in^ 

(E minimized) is given by 


jgF(t) = ^ A n L n (st); 


n= 1 


where 


A n (s) z | F(t)L n (st)dt 


( 20 ) 


( 21 ) 


The integral squar 


re error is now given by 


F(t)^dt 


2,n. - 


V- 2 

> An > 0 

nrrl 


( 22 ) 


E is minimized by choosing s such that ^ A n 1 


is maximum. 


n= 


Theorem 7, completeness theorem: 

If 


26 
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"SV! w^st.rr rt. vM«i; 


OO 

J F(t) 2 dt 

0 

is finite, and the Laplace transform of F(t) f f(s) exists, then 
E + 0 as N + 03 


Theorem 8: 

Let the Laplace transform of F(t) be 


f(s) = J F(t)e- St dt 
0 


(23) 


Then 


n 

A r ,(s) = 2 ^ n a i 1 ( is ) 

i=1 


(2*0 


Theorem 9: 


^F(t) can be written as 


tjF(t) = 


N 

y NB n *“ nat 

n- 1 


(25) 


.where 


N'-’n 


/•. r v ^ \ , v.. r / 7 ~ \ 

<-* N l -’n1 u :- y r T hh'nj 1 * 


4 >rV> . 

* (v-rnv 4 v 4 "-" | 


?r, \ 


27 





and where 


H b ij s N b ji = — 'yT' k a i k a 1 (k a m = 0 for m > 
2s ' 


k) 


k= 1 


or 
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(27) 


1 . (-1) 1+ J 1 

N b ij = N b ji 


1 (N + i)l (N + j)! 


2(i -*• j) i!(i - 1)! j ! ( j - 1)! (N - i)J (N - j)! 


(28) 


min a 2: 


Y N b ij / (-1) N ~J (N + j)< 

(n~ j ) s j ! 77 “^*~TTT 


i: 1 


x + i 2 


(x - 1)(x - 2)- • • (x - N) 1 

(x ■*• i ) ( x + 2) - • • (77 n7 —■ 


(29) 


Theorem 10: 


. 1 v N °n 

f( is) = - \ 

s i + n 

ns 1 


(30) 


Theorem 1 1 : 


y £. 


n=1 


N 

N N 13 n f( ns) 
ii" 1 


(3D 


where ^jB n was given by equation .26. 


28 



1 


Theorem 12: 
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^F ( 0 ) = \/2ns A n 

£-4 f 

n=1 


(32) 


Theorem 13: 


G(t)L ra (st)L n (st)dt ~j ^ m a j n8 


SJ 


U s i 


where 


n 


n^j 


5 j = > n a ifif(i + J)s 


.. .-'f 1 

x^ “ A l 
i - 1 


(33) 


(3*0 


where g(s) is the Laplace transform of G(t). 

Theorem 14: 

The best approximation to the jth derivative of F(t) is 

■ N" 

N FU)(t.) = Y jA n L n (st) _ ' (35) 

n=1 


where 


29 ’ 
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jA n = ]> n a 1 <la)JrCls) - F(+0) £ n aiUs)l-1 


isl 


1=1 


(36) 


dF 

dt 


^ n a i ( is )^ 2 - 


t=+0 i = 1 


d j" 1 F 


d 2 F 

dt 2 


]T n ai(i a)J-3 


t=+0 


L=1 


dtJ-1 


n 

I -* 


t=+0 


i= 1 


Note 


HP ( J , .(tj 


dJ K F(t) 
d t J 


(37) 


For example, IF j = 1, the first derivative, then 


1 A n = N n aj_( is) f(is) » /?sn F( -#-0 ) 
i = 1 


(38) 


Note equation 10, corollary B, 


n s 
i = 1 


was used to obtain equation 38. The value of F(-t-O) can be obtained f rom 
the initial value theorem. 


H>U) = Lira *f{s) 

S ■+ <■» 


( 39 ) 


30 







If j s 2 t the second derivative, then 
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? A n - y n a^( is)2f( is) - /Ssn n* : sF(+0) -/23a — 
«£-* * dt 

i-1 


t=+0 


(^0) 


If j = 3 


« •*. x\ d 

3 A n - / n a i^ is ^f^ s ^ - /2sn — (n^ + 1)s^F(+0) 

*— * 2 

1 = 1 


., dF 

/2sil n'-c 

dt 


- /?. sn 


d 2 F 

dt* 


t=+0 


t =+0 


('41) 


For j = ^ 


5. ,m * n c 

i|A n = \ n a ii>)^F( is) - /2sn —(n^ + + 1 )s3f(+0) 

L~t 6 

i = 1 


1 dF 
- /?IFT n2(n? + 1) S 2 — 

2 dT 


nr — * p 

/ 2 s n n"*s 


d^r 


dt" 


t-=+0 


t =+0 


/ 2sn 


d"F 

dt3 


t~4 0 


W) 


Theorem 15: 


If 



(M> 



. 80FK20 


ther 


A n = (~l) n +U /5ns 


(s - a) (2s - a) f ** ((n - 1 )s - a) 
(s + a) (2s + a)*** (ns + a) 


<M) 


A ] s A /5s" — 


s + a 


(*5) 


Note the results for A s 1 and a = 0, F(t) a unit step function.. In this 


case 


A n = C-1)n +1 


/5nT 


Hence 


(t6) 


n fu) 


v \ (-D n+1 = 

\ — L n ( st ) 
. /2 ns 
nsl 


W> 


From corollary D, L n (0) = /5hs t hence 

N even (^8) 

N odd 

The equations shown in theorem 15 'are useful for testing the accuracy of 
computer computations. 


nF(0) = 0 

r 2 




